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Abstract 

We consider the description of the fractal media that uses the fractional integrals. 
We derive the fractional generalizations of the equation that defines the medium mass. 
We prove that the fractional integrals can be used to describe the media with noninteger 
mass dimensions. The fractional equation of continuity is considered. 
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1 Introduction 

The application of fractals in physics IH| is far ranging, from the dynamics of fluids 
in porous media to resistivity networks in electronics. The cornerstone of fractals is the 
meaning of dimension, specifically the fractal dimension. Fractal dimension can be best 
calculated by box counting method which means drawing a box of size R and counting the 
mass inside. Fractal models of media are enjoying considerable popularity. This is due in 
part to the relatively small number of parameters that can define a fractal medium of great 
complexity and rich structure. 

Derivatives and integrals of fractional order have found many applications in recent stud- 
ies of scaling phenomena [3 IU |H1 E| • It is interesting to use fractional integration 
to consider the properties of the fractional media. It is interesting to find the connection 
between the fractional integrals and fractals. The natural questions arise: what could be 
the physical meaning of the fractional integration? This physical meaning can be following: 
the fractional integration can be considered as an integration in some fractional space. In 
order to use this interpretation we must define a fractional space. The first interpretation 
of the fractional space is connected with fractional dimension space. The fractional dimen- 
sion interpretation follows from the formulas for dimensional regularizations. If we use the 
well-known formulas for dimensional regularizations [TT], then we get that the fractional 
integration can be considered as an integration in the fractional dimension space up to 
the numerical factor T(D /2) / ' (2ii D ^ 2 T(D)). This interpretation was suggested in Ref. [TU] . 
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In this Letter we consider the second physical interpretation of the coordinate fractional 
integration. This interpretation follows from the fractional measure of space (TOl that is used 
in the fractional integrals. We consider the mass fractional dimension and the fractional 
generalizations of the equation that defines the mass of the medium. We prove that the 
fractional integrals can be used to describe fractal media with noninteger mass dimensions. 
The fractional integrals can be used not only to calculate the mass dimensions of fractal 
media. Fractional integration can be used to describe the dynamical processes in the fractal 
media. Using fractional integrals, we can derive the fractional generalization of the dynamical 
equations. The fractional generalization of the Liouville equation was suggested in Ref. 
1CL Using the Fourier transform, we can introduce fractional dynamical equations with 
coordinate fractional derivatives. 

In Section 2, the definition of mass dimension and the fractional generalization of the 
medium mass equation is considered. In Section 3, we consider the properties of the fractal 
media. We define fractality and homogeneity properties of the media. In Section 4, we 
discuss the local density of the fractal media. In Section 5, the fractional generalizations 
of the equations are considered. We prove that these equations can be used to describe 
fractal media. In section 6, we consider the fractional equation of continuity. Finally, a 
short conclusion is given in Section 7. 

2 Mass Fractal Dimension 

Equations that define the fractal dimensions have the passage to the limit. This passage 
makes difficult the practical application to the real fractal media. The other dimensions, 
which can be easy calculated from the experimental data, are used in empirical investigations. 
For example, the mass fractional dimension [T3*l ITi] can be easy measured. 
The properties of the fractal media like mass obeys a power law relation 

M(R) = kR Dm , (1) 

where M is the mass of fractal medium, R is a box size (or a sphere radius), and D m is a 
mass fractal dimension. Amount of mass of a medium inside a box of size R has a power 
law relation (fTJ) . 

Fractal dimension can be calculated by box counting method which means drawing a box 
of size R and counting the mass inside. To calculate the mass fractal dimension, take the 
logarithm (JTJ) both sides 

ln(M) = D m ln(R) + Ink. 

Log-log plot of M and R gives us the slope D m , the fractal dimension. When we graph 
ln(M) as a function of ln(R), we get a value of about D m which is the fractal dimension of 
fractal media. 

The power law relation (0) can be naturally derived by using the fractional integral. 

In order to describe the fractal media, we suggest to use the space with fractional measure. 
In this case, we can use the constant density distribution for homogeneous fractal media. If 
we use p(x) = po = const, then we get (JH). In the next sections, we prove that the mass 
fractal dimension is connected with the order of fractional integrals. Therefore the fractional 
integrals can be used to describe fractal media with non-integer mass dimensions. 
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3 Fractal Media 



Let us consider the region Wa in 3-dimensional Euclidean space E 3 , where A is the midpoint 
of this region. The volume of the region Wa is denoted by V(Wa). If the region Wa is 
a ball with the radius Ra, then the midpoint A is a center of the ball, and the volume 

V(W A ) = (4/3)7^ . 

The mass of the region Wa in the fractal media is denoted by M(Wa)- Let p(W A ) is an 
average mass density of the region Wa of the fractal medium. This density is defined by the 
equation 

p(W A ) = M(W A )/V(W A ). 

In the general case, the fractal media cannot be considered as continuous media. There 
are points and domains that do not filled by the medium particles. These domains can be 
called the porous. The fractal media can be considered as continuous media for the scales 
much more than mean value of the pore size R p . 

The fractal media can be characterized by the following property: 

// the volume V(Wa) of the region Wa increase, then the average mass density p(Wa) is 
decrease. 

This property is satisfied for all points A and all regions Wa if the volume of region is 
much more the average value of porous volume (V(W,a) >> V p ). For the ball region of the 
fractal media, this property can be described by the power law 

p(W A ) = po(R/R p ) d - 3 , 

where p is a constant value; R is the radius of the ball Wa- Here R >> R p , where R p is a 
mean radius of the porous sphere. 

The fractality of medium means that the mass of this medium in any region of Euclidian 
space E 3 increase more slowly than the volume of this region. For the ball region W, we 
have 

M(R) ~ R D , D < 3. 
The characteristic property of fractal media can be considered in the following forms: 

For all regions Wa and Wb in the fractal media such that Wa C Wb 
and V(Wa) < V(Wb), we have that the corresponding average mass densities satisfy the 
inequality p(Wa) > p(Wb), i.e., 

W A CW B , V(W A )<V(W B )=>p(W A )>p(W B ). (2) 

We would like to consider the fractal medium that have some homogeneous property. 
Now we shall give the definition of homogeneous fractal media. 

The fractal media are called homogeneous if the following property is satisfied: 

For all regions Wa and Wb of the homogeneous fractal media such that the volumes are 
equal V(Wa) = V(Wb), we have average densities of these regions are equal too 
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p(W A )=p(W B ), i-e-, 



V(W A )=V(W B ) => p(W A )=p(W B ) 



(3) 



Fractal medium is called a homogeneous fractal medium if the average density value of 
the region does not depends on the translation of this region. 

Note that the wide class of the fractal media satisfies the homogeneous property. In 
many cases, we can consider the porous media [T3J HE] , polymers [T7] , colloid agregates [TH] , 
and aerogels JH] as homogeneous fractal media. The dendrites cannot be considered as a 
homogeneous fractal medium. 

We can generalize the fractal property (J2J) for the homogeneous fractal media. This gen- 
eralization is connected with consideration of two different points and two different regions 
of the media. For the homogeneous fractal media we can remove the restriction Wa C Wb 
in Eq. (J2J). These regions can satisfy the condition Wa H Wb = 0- Therefore we have the 
following property: 

For all regions Wa and Wb of the homogeneous fractal medium in the Euclidean space E 3 
that satisfies the inequality V(Wa) < V(W B ), we have that the average densities are con- 
nected by the inequality p(Wa) > p(W B )- 

As the result, we have the following properties for all regions Wa and W B in the Eu- 
clidean space E 3 with homogeneous fractal media: 



(1) Fractality: if V(W A ) < V(W B ), then p(W A ) > p(W B ). 

(2) Homogeneity: if V(W A ) = V(W B ), then p(W A ) = p(W B ) and M{W A ) = M(W B ). 

4 Local Density of Fractal Media 



In the general case, the fractal media cannot be considered as continuous media. There are 
points and domains that do not filled of the medium particles. These domains can be called 
the porous. The fractal media can be considered as continuous media for the scales much 
more than mean value of the pore size R p . In this case, we can use the integration and 
differentiation. 

In many cases, the fractal media are described by the equations with the integer integra- 
tion. However, this description can be incorrect for fractal media. For example, the mass of 
the medium is derived by the following equation 



If we describe local density p(r) as the single-valued function, then we cannot use Eq. (j3J) 
with the integer integral for the homogeneous fractal media. 

In order to prove this statement, we consider the points A, B, and C, such that Ra = 
\AC\ > Rb = \BC\. Obviously, we require that Ra » R p and R B » R p . Let us describe 




(4) 
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the mass of the media using the integer integrals. The mass of the regions Wa and Wb of 
the homogeneous fractal media are defined by the equations 

M(W A ) = ! PA(r)rf 3 r, M(W B ) = f p B (r)d 3 r. (5) 

JWa JW b 

We use the index A and B in the local densities, since we would like to prove that Pa( v c) 7^ 
Ps( r c*)- Obviously, that the local densities pa and pB of the continuous medium satisfy the 
condition p A (r c ) = p B (r c ). 

In order to satisfy the fractal property of the medium and Eqs. (JSJ), we must use the 
following local densities: 

PaOO = Po(|r - r A \/R^ D 3 , p B {r) = p (|r - t b \/R^ D 3 . (6) 

Let us consider the point C such that \AC\ > \BC\. Using Eq. (JHJ), we have the local 
density in this point in the form 



PA(r c ) = Pv(\AC\/R p ) D ~\ p B (r c ) = p Q (\BC\/ R p )°^ . 



Since \AC\ > \BC\, it follows that Pa(?c) < Pb(?c)- Obviously, that we have two different 
values for one point C of fractal medium. Therefore we cannot use Eq. (J3J) with the integer 
integrals for the homogeneous fractal media. We can use Eq. (0J) with the local density 
p(r) ~ |r| D_3 only for the non- homogeneous distribution of the medium mass that has the 
single out point. This point is defined by r = 0. 

As the result, we have that the fractality and homogeneity properties of the medium 
cannot be satisfied simultaneously by Eq. (|IJ) with the integer integrals. 

In order to satisfy these properties, we must use the generalization of Eq. (jlj such that 
fractality and homogeneity properties can be realized in the form: 

(1) Fractality: the mass of the region W of fractal medium obeys a power law relation 

M W ^M (2M!) D , (7, 

JXp 

where D < 3. For the ball region Wa, we have 

M(R) = M Q (R/R P ) D , 

where R — \ r — r a \ ■ 

(2) Homogeneity: the local density of homogeneous fractal medium is translation invariant 
value that have the form 

p(r) = po = const. (8) 



We can realize these requirements by the equation 

M(W) = [ p(r)dp D (r), 

JW 
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where /xp is a new measure of the space E 3 . Using Eq. (JBJ), we have 



M(W) = po [ dfi D (r) = p V D (W), 



where Vb(W) is the volume of the region W. Therefore we get that the homogeneous fractal 
media are described by the measure space (E, p) such that 



where V(VK) is the usual volume of the region W, i.e., we have Vd(W) ~ (V{W)) D ^ 3 . If we 
consider the ball region, we get 



where D < 3, and Vq = (47r /3) d ^ 3 (po/Mq). In the next section, we prove that the natural 
generalization of Eq. (jlj) uses the fractional integration. 

5 Fractional Equations 

Let us prove that the fractional integration allows us to generalize Eq. (@J) such that the 
fractality and homogeneity conditions are realized in the form 



Let us define the fractional integral in Euclidean space E 3 in the Riesz form ^2] by the 




V D (W) = V (R/R P ) 



M(W) ~ (V(W)) D/3 , V D (W) ~ (V(W)) D/3 , p(r) = const. 



equation 




(9) 



where r G W, and we use 




2 D n 3 / 2 T(D/2) 
T(l/2) 



Using R = r — r , we can be rewritten Eq. Q in an equivalent form 



(I D p) (r ) =73^ (D) I p(R + r )|Rr- 3 d 3 R. 
Jw 



Note that 73(3 — 0) ^ 1, where 



7s(3 - 0) 




2 3 7r 3 / 2 r(3/2) 

r(i/2) 



Using notations (0), we can write Eq. (jlj) in the form 



M(W0=7 3 (3-0)(/ 3 p)(r ). 
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The fractional generalization has the form 

M(^)= 73 (3-0)(/ z V)(r ). 

We use the factor 7 3 (3 — 0) to derive Eq. (jlj) in the limit D — > 3 — 0. The fractional equation 
for M(W) can be written in the equivalent form 

> 3 - D r(3/2) 

V(D/2) Jw 

where we use 

73(3-0) 73 -'( O ) . 

Here and late we use the initial points in the integrals are set to zero. 
If we have p(r) = p = const and the ball region W, then 



M(W) = / w p(R + r Q )|Rr- 3 d 3 R, (10) 

2 3 - D r(3/2) 



M(W) = p 23 '? r(3 ( 2) / |R| D - 3 rf 3 R- 

Using the spherical coordinates, we get 

_ ^ 5 ^r(3/2) r _ 2^^(3/2) 

M(W) " T(D/2) P ° Jw R dR ~ DT(D/2) P ° R 



D 



where R = |R|. As the result, we have M(W) ~ R D , i.e., we derive Eq. ((7|) up to the 
numerical factor. 

We can suppose that Eq. (fTUj) can be used for non-homogeneous fractal media. In this 
case, the distribution depends on R = |R|, and the angels. 

Let us consider the spherical- homogeneous fractal media that have the single out point. 
These media are defined by the following properties: 

(1) For all regions Wa and Wb in the Euclidian space E 3 such that \OA\ = \OB\ = R and 
R much more than the size of the regions (R 3 « V(Wa) = V(Wb) « R 3 ), we have the 
relation p(Wa) = p{W B )- 

(2) For all regions Wa such that the size of the region is much smaller than distance \OA\ = R 
(Rp << V(Wa) « R?), we have that the density p(Wa) depends on the distance R = \OA\. 

For example, we can consider the power law in the form 

p(W A ) ~ \OAf. (11) 

The example of the spherical-homogeneous fractal media we can be realized in nature as 
dendrites. 

If we consider the non-homogeneous fractal media, then the local density p = p(R) have 
the single out point. For the spherical symmetric case, we have p(R) = p(|R|) and the 
equation can be represented in the form 

vr2 5 - D r(3/2) fR 
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For the simple case we have the distribution in the form 



p(R)=p(|R|) = c |R| /3 . 
Therefore the mass of fractal medium in the region W is 

M (W) = c 23 "" >r ( 3 /2) r | R |D +/3 -3 d 3 R _ 

V ; r(D/2) Jw l 1 

Using the spherical coordinates, we have 

M(W) - C0 T(D/2) L R dR ~ (D + f3)r(D/2) R ' 

In this case, the mass- dimension is equal to D m = D + (3. 

To calculate (3, we can use the definition of the spherical-homogeneous fractal media and 
the generalization of the box counting method. This generalization must use the box regions 
W that satisfy the conditions ^(W 7 ) << R from the definition of spherical-homogeneous 
fractal media. 



6 Fractional Equation of Continuity 

The fractional integrals can be used not only to calculate the mass dimensions of fractal 
media. Fractional integration can be used to describe the dynamical processes in the fractal 
media. Using fractional integrals, we can derive the fractional generalization of dynamical 
equations [HI. Let us derive the fractional analog of the equation of continuity for the 
fractional media. 

Let us consider a domain Wq for the time t = 0. In the Hamilton picture we have 



/ p(Rt,t)dii D (R t ) = [ p(Ro, O)d^(Ro) 

JWt JWn 



'Wt JW 

where we use the following notation 

^ D (R) = 2 3 ^r(3/2) D 
h- UK ) T(D/2) 1 1 

Using the replacement of variables R t = R 4 (R ), where R is a Lagrangian variable, we get 
/ p(R t ,t)\R t \ D - 3 ^dR = [ p(R ,0)|R | D - 3 ciRo, 

where dRj/SRo is Jacobian. Since Wq is an arbitrary domain we have 

p(R t ,t)dfj, D (Rt) = p(Ro,0)rf/i D (R ), 

or the equivalent form 

p(R t ,t)\R t \ D - 3 ^ = p(R ,0)|R |r 3 . 
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Differentiating this equation in time t, we obtain 

<jp(R t ,t) 
eft 

where rf/rft is a total time derivative 



n D (Ri,t)p(Ri,t) =0, (12) 



d d dR t d 
+ 



dt dt dt dR t 

The function 

n D (R t ,t) = ^-ln(\R t \ D - 3 ^ 
dt v oR ' 

describes the velocity of volume change. Eq. (|T2*|) is a fractional Liouville equation in the 
Hamilton picture. Using the equation 



dR t 
~dT 

we get the function 



V*(R,t), 



As the result we have 



fiD(R ' t} - IBP + SR' 



where we use <9|R t |/<9R t = R t /|R t |. 

In the general case (D ^ 3), the function fl D is not equal to zero (Qd ^ 0) for solenoidal 
velocity fields (divV = 0). If D = 3, we have Qd ^ only for non-solenoidal velocity field. 
Therefore the flow in the fractional medium is similar to the non-solenoidal flow. 

As the result, we have the following equation of continuity for fractal media: 

dp(R t ,t) d P (R t ,t) (D-3)(Rt,V f ) ^ V '-n 

~^r~ + ~dkv~ + — \Rj* — p{ } + p( } 9R t ~ °- 

For the homogemeous fractal media, we have p = const and the equation of continuity leads 
us to the equation 

(D-3)(R t ,V f ) gVe 

IRI 2 9R t • 

Therefore, we get the non-solenoidal flow of the velocity divV = d\/dR ^ 0. 
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7 Conclusion 



The application of fractals in physics is far ranging, from the dynamics of fluids in porous 
media to resistivity networks in electronics. The cornerstone of fractals is the meaning of 
dimension, specifically the fractal dimension. Fractal dimension can be best calculated by 
box counting method which means drawing a box of size R and counting the mass inside. 
When we graph ln(M) as a function of ln(R) we get a value of about D m which is the 
fractal dimension of fractal media. The mass fractal dimensions of the media can be easy 
measured by experiments. The experimental determination of the mass fractal dimension 
can be realized by the usual box counting methods. The mass dimensions were measured for 
porous materials [13113, polymers [T7], and colloid agregates [TH]. For the porous materials 
the fractal dimension is proportional to porosity of the medium. Fractal models of porous 
media are enjoying considerable popularity [201 |2U 123 E3 EH I2S]- This is due in part to 
the relatively small number of parameters that can define a fractal porous medium of great 
complexity and rich structure. 

We suppose that the concept of fractional integration provides an alternative approach 
to describe the fractal media. In this Letter we prove that mass fractal media can be 
described by using the fractional integrals. We consider the fractional generalizations of the 
equation that defines the medium mass. This fractional generalization allows us to realize 
the consistent description of the fractal mass dimension of the media. 

The fractional integrals can be used in order to describe the dynamical processes in the 
fractal media. The physics of fractal media can be successful described by the fractional 
integration. Fractional integration approach is potentially more useful for physics of fractal 
media than traditional methods that use the integer integration. 

Note that the Liouville equation is a cornerstone of the statistical mechanics [213 EH EH] • 
The fractional generalization of the Liouville equation was suggested in Ref . ^U] • This frac- 
tional Liouville equation was derived from the fractional generalization of the normalization 
condition JU] that uses fractional integrals. The fractional generalization of Liouville equa- 
tion allows us to derive the fractional analogs of Bogoliubov equations [29, , Fokker-Planck 
equation, Enskog transport equation, and hydrodynamic equations. Using the Fourier trans- 
form, we can get the fractional dynamical equations with coordinate fractional derivatives. 
Using the methods suggested in Refs. |30J EU E2] , we can realize the Weyl quantization of 
the suggested fractional equations. 
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